Dynamical Decoupling of Qubits in Spin Bath under Periodic Quantum Control 
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We investigate the feasibility for the preservation of coherence and entanglement of one and two 
spin qubits coupled to an interacting quantum spin-1/2 chain within the dynamical decoupling 
(DD) scheme. The performance is examined by counting number of computing pulses that can be 
applied periodically with period of T before qubits become decoherent, while identical decoupling 
pulse sequence is applied within each cycle. By considering pulses with mixed directions and finite 
width controlled by magnetic fields, it is shown that pulse-width accumulation degrades the perfor- 
mance of sequences with larger number of pulses and feasible magnetic fields in practice restrict the 
consideration to sequences with number of decoupling pulses being less than 10 within each cycle. 
Furthermore, within each cycle T, exact nontrivial pulse sequences are found for the first time to 
suppress the qubit-bath coupling to 0(T N+1 ) progressively with minimum number of pulses being 
4,7,12 for N = 1,2,3. These sequences, when applied to all qubits, are shown to preserve both 
the entanglement and coherence. Based on time-dependent density matrix renormalization, our 
numerical results show that for modest magnetic fields (10-40 Tesla) available in laboratories, the 
overall performance is optimized when number of pulses in each cycle is 4 or 7 with pulse directions 
be alternating between x and z. Our results provide useful guides for the preservation of coherence 
and entanglement of spin qubits in solid state. 

PACS numbers: 03.65.Yz, 03.65. Ud, 03.67.-a 



I. INTRODUCTION 

The dream of building quantum computers has driven 
intensive investigations on quantum information process- 
ing during the past decade. Nonetheless, due to the 
ubiquitous decoherence problem, the progress made so 
far has been limited. Since the processing of a real 
quantum system includes inevitable disturbance from 
the outside world, the central challenge is to find ways 
to control or even eliminate the decoherence. There 
are different strategies proposed to overcome decoher- 
ence such as dynamical decoupling (DD) 1-7 , quantum 
error correction 8 - - — , and decoherence-free subspaceii~— . 
While different strategies have their own advantages, the 
dynamical decoupling represents the oldest effort along 
this direction and have been known as a mature tech- 
nique employed in Nuclear Magnetic Resonance (NMR) 
experiments. Theoretically, it has been rigorously shown 
that DD provides upper bounds for error of reduced den- 
sity matrix caused by quantum evolution^. Recent NMR 
experiments further indicate that dynamical decoupling 
does preserve coherence of a nuclear-spin qubit 15 . These 
facts clearly indicate that DD is promising in providing 
a practical solution to defeating decoherence. 

To implement the scheme of dynamical decoupling, ex- 
plicit pulse sequence has to be constructed. Various pulse 
sequences were proposed and developed. Hanh's spin 
echo (SE)i& and Carr-Purcell-Meiboom-Gill (CPMG)i2 
were brought up in the beginning. Later, concatenated 
dynamical decoupling (CDD)S sequence and Uhrig's dy- 
namical decoupling (UDD) 7 sequence were proposed. 
With so many pulse sequences available, one still needs 
to address the central issue in the scheme of dynamical 



decoupling: what is the sequence that has the best perfor- 
mance in suppressing decoherence while viable quantum 
manipulations are kept? The issue has been addressed 
by considering a given cycle of T in which pulses are ap- 
plied. The performance is examined by number of pulses 
needed for suppressing the qubit-bath coupling to the 
order 0(T m ). When durations of pulses are ignored, it 
was recently shown by Yang and Lir*i£ that for a single 
qubit interacting with bath with Ising-like coupling, the 
UDD-N pulse sequence can suppress the pure dephasing 
to 0{T N+1 ). However in addition to the control of de- 
phasing, one also needs to control longitudinal relaxation. 
This would be necessary when the coupling between the 
qubit and the spin bath is Heisenberg-like. In this case, 
Yang and Liu 18 showed that the UDD-N pulse sequence 
can not eliminate the longitudinal relaxation and the de- 
phasing to 0(T N+1 ) at the same time. This calls for a 
closer examination on the minimum number of pulses for 
suppressing the qubit-bath coupling to the desired order 
0(T N+1 ). Recently, a quadratic DD sequence (QDD) 
that concatenating x-direction and z-direction UDD se- 
quences is proposed™ Although QDD is shown to sup- 
press general decoherence to 0(T N+1 ) by using (N + 1) 2 
pulse intervals, the sequence is not optimal and the issue 
of finding the optimal sequence still remains. 

From theoretical point of view, if the qubit-bath cou- 
pling is the only Hamiltonian that governs qubits, the 
reduced density matrix p r {T) of qubits includes all un- 
desired dynamics. Therefore, for a given T and order 
N, a sequence is optimal if it suppresses all operators in 
p r {T) to 0(T N+1 ). However, in order to perform com- 
puting, one also needs a strategy for inserting comput- 
ing pulses 2 ^. It is clear that if the suppression due to 
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DD pulses is indiscriminating, the desired dynamics due 
to computing will also be suppressed unless computing 
pulses form another commuting DD-pulse sequence. In 
this case, one decomposes DD pulses into cycles separat- 
ing by computing pulses. The performance of DD pulses 
is then examined by number of computing pulses that 
can be applied before the system becomes decoherent. In 
addition to the issue of how to insert computing pulses, 
the finite duration of pulses also represents an important 
constraint. The accumulation of pulse- width is seen to 
degrade the performance of DD pulses 15 . It is therefore 
important to compare performance of sequences with dif- 
ferent orders. So far, most construction of DD pulses fo- 
cuses on single qubit. It is known that the entanglement 
is particularly important for characterizing the quantum 
state of multi-qubits and plays the crucial role in quan- 
tum information processing. There have been a few inves- 
tigations of effects of DD pulses on multiqubits. West et 
al. investigated fidelity of quantum states of four nuclear 
spin-qubits in the decoherence free spaced and found DD 
does preserve the fidelity. There have also been studies 
based on pulse control of the entanglement for two qubits 
in rather simplified models^—. Nonetheless, it is still 
not clear what would be an optimized sequence for pre- 
serving the entanglement. 

In this paper, we investigate the feasibility for preser- 
vation of decoherence and entanglement of spin qubits 
within the DD scheme in solid state system. One or two 
spin qubits coupled to an interacting quantum spin- 1/2 
chain are considered. We shall examine different strate- 
gies for inserting computing pulses within each cycle and 
demonstrate that computing after decoupling performs 
the best. We then examine the feasibility by insert- 
ing computing pulses periodically with period T within 
which the same dynamical decoupling pulses are applied. 
It is shown that error induced by pulse-width accumula- 
tion restricts the consideration to sequences with num- 
ber of pulses being less than 10 within each cycle. Fur- 
thermore, within each cycle T, exact nontrivial pulse se- 
quences can be constructed to suppress the qubit-bath 
coupling to 0{T N+1 ) progressively with number of pulses 
being 4,7,12 for N = 1,2,3. Based on time-dependent 
density matrix renormalization (t-DMRG), our numeri- 
cal results show that for modest magnetic fields (10-40 
Tesla) available in laboratories, the overall performance 
is optimized when number of pulses in each cycle is 4 or 
7 with pulse directions be alternating between x and z. 

This paper is organized as follows. In Sec. II, we 
present our model Hamiltonian and briefly discuss how 
to apply t-DMRG to analyze the model Hamiltonian. 
We shall outline the general framework for calculating 
the dephasing and longitudinal relaxation. In particu- 
lar, we point out that for general coupling between the 
qubit and the bath, the preservation of either coherence 
or entanglement is determined directly by the evolution 
operator U(t). In Sec. Ill, we analyze decoherence and 
longitudinal relaxation of a single qubit by considering 
pulses with mixed directions. We will explicitly construct 



pulse sequences for suppressing lower orders of U (T) up 
to 0(T 4 ). Furthermore, different strategies for inserting 
computing pulses are compared. We find that comput- 
ing after decoupling performs the best. Therefore, we 
extend DD over a cycle of T to the periodic scheme in 
which computing pulses are inserted at nT. By consid- 
ering the finite duration of pulses, we further analyze 
dynamics defined at nT. We show that for available mag- 
netic fields, the number of quantum manipulations can be 
maximized by using sequence consisting of 4 or 7 pulses. 
Sec. IV is devoted to investigate the entanglement of 
two qubits. We shall show that regardless whether two 
qubits are strongly entangled or non-entangled, the en- 
tanglement can be preserved by using the same sequence 
that suppresses the decoherence of a single qubit. We 
further show that for general multi-qubits scenario, en- 
tanglement can be preserved by applying the same se- 
quence to all qubits if separations between qubits are 
sufficiently large. In Sec. V, we summarize our results 
and discuss possible generalization to dynamically decou- 
ple multi-qubits from the environment. In Appendix A, 
we explicitly construct equivalent sequences for N — 2. 



II. THEORETICAL FORMULATION AND 
GENERAL CONSIDERATION 

We consider a system-bath model which is described by 
the total Hamiltonian Hq = H sys + H{, a th + Hint , where 
H sys is the Hamiltonian of a single or two qubits system, 
Hbath is the Hamiltonian of a spin bath and Hi nt repre- 
sents the interaction between qubits and the bath. The 
system Hamiltonian H sys is generally zero unless com- 
puting pulses or decoupling pulses are applied. Generally, 
computing pulses can be also spread over all times^. In 
this case, one has 

H S y S ^s-B, (1) 

where s is the qubit spin operator and B is the corre- 
sponding magnetic field for computing. The spin bath is 
a spin chain generally characterized by the XX Z Heisen- 
berg model 

'Hbath = J ^ {SfSf +1 + + A-SfSf+i) , (2) 

where J > 0. It is known that the XXZ Heisenberg 
model has a very rich structured The decoherence and 
entanglement dynamics induced by such kind of spin bath 
have been recently investigated ! 26 ' 27 In order not to be 
masked by dynamics of the order parameter—, we shall 
focus on the XY regime where |A| < 1. Since the case 
of A / behaves qualitatively the same as A = case, 
we shall simply set A = with understanding that our 
results are also applicable to A / 0. Two specific forms 
of the qubit-bath coupling Hamiltonian Hi nt are consid- 
ered. In the first scenario we consider Ising-like coupling 
Hi„t = cs z S* for a single qubit or H int — e(sfSf + s^S*) 
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for two qubits, where s z (Sf) is the qubit (spin bath) op- 
erators and is the single site of spin chain to which 
the qubit is coupled to. It gives rise to pure dephasing 
of the qubit. To minimize the boundary effect in numer- 
ical calculation, is usually taken to be the central 
site of the chain. In the second scenario we consider 
Heisenberg-like coupling Hi nt = es- Si for a single qubit 
or Hi n t = e(si • Si + S2 ■ Sj) for two qubits, which gives 
rise to both dephasing and energy relaxation. 
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FIG. 1: Illustration of a dynamical decoupling pulse sequence 
in the periodic decoupling scheme. Here a sequence of pulses, 
each with width t p centered at tk, are applied periodically 
with period T. The directions of pulses are alternating and 
to compensate the even/odd effect of pulse direction, a par- 
ity pulse marked by slashed line is applied at the end of the 
sequence. At t* = nT, computing pulses of width r for 
quantum processing are applied. For convenience, we denote 
intervals between centers of pulses by anT. 



H c (t) 



i-k 



A(t - t k )si 



(5) 



Here Si represents the spin of i th qubit and h], is the 
control axis during the fc th pulse that is applied at the 
« th qubit. A(t — tk) is a square function centered at tk 
with width t p . The magnitude of A is HbB with being 
the Bohr magneton and B be the magnetic field so that 
Jt p = 107T J(meV)/_B(Tesla). Experimentally, accessible 
magnetic field B will impose a lower bound on the pulse 
duration r p . 

We shall focus on periodic dynamical decoupling where 
H c (t + T) = H c {t) and T is the period of the pulse se- 
quence. Manipulating or computing pulses with width 
r are applied at t* = T, 2T, 3T, • • • . Since both t q and 
t p are determined by available magnetic fields, we shall 
assume t q = t p . Furthermore, in order that pulses are 
non-overlapping, one requires 



T > max 



(6) 



Since o^min) < 1/(N + 1), we have T > (N + l)r p . 
Therefore, T depends on JV. In the following we will 
denote T by T/v to indicate its explicit dependence on 
N. Formally the evolution of the total system is dictated 
by the evolution operator 



(7) 



For most of our numerical work, we shall focus on the 
initial state in which the total state of the system is a 
product state of the form: 



|*(0)> = \g) ® \G), 



(3) 



where qubits are in some particular state of interest \g) 
while the bath is in its ground state \G). Our results, 
however, are based on the consideration of evolution op- 
erator directly. Hence similar results are also found for 
other initial states. 

A pulse sequence in a cycle of T may contain N de- 
coupling pulses centered at tk where k = 1, 2, • • ■ , N as 
illustrated in Fig. [TJ Because number of pulses along 
the same direction may not be even, to compensate this 
parity effect, a parity pulse is added at the end of the 
sequence. We shall denoted intervals between tk by a{T. 
Typically in an experiment, these pulses are generated 
by a magnetic field B{t)h and have the same width t p 
. We concentrate on the pulse sequences in which each 
pulse gives rise to a n rotation along a given control axis 
n for some qubit. If pulses are characterized by a time- 
dependent control Hamiltonian H c (t), the total system 
is then characterized by the Hamiltonian 



H(t) = H c (t) + H , 



(4) 



where we have set h to one and T is the time-ordering op- 
erator. During each pulse period, the evolution operator 
for the kth pulse can be written as 



Te 



- i: k k -, p/ 2 Ho{s)ds -i /^_+;;/ 2 2 H c (s)ds 



xT e 



t k +r p /2 



H (s)ds e O(NejT^) 



(8) 



where the second equality introduces an error of the order 
Tp. Since for each qubit one has J^^ P /2 H c ds — nsi ■ h k , 
the evolution operator can be expressed as 

t/W(t) = {-i) N e - lHo{t - tN) ii t d l ■ fe N e- lHa{tN - tN - i) x 

[1 + 0{NejT 2 p )\ , (9) 



where a = 2s are the Pauli matrices for the qubit and the 
rotational matrix e -' l ^i s - n k) [ s reduced to ia ■ hk- Eq.© 
thus implies that to the order of T p , N pulses of finite 
width can be consider as ideal pulses without width so 
that qubits arc flipped right after tk- 

In general it is difficult to exactly evaluate U N (t). t- 
DMRG, however, provides a way to efficiently evolve such 
a state with high accuracy for a quasi-one dimensional 
system. We first use static DMRG to find the ground 
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state of the spin chain bath for a given A and then use 
the method of t-DMRG to evaluate U N (t)\${0)) numer- 
ically. We note that the degrees of freedom of the qubits 
are kept exactly during the t-DMRG calculation by tar- 
geting an appropriate state. The dimension of the trun- 
cated Hilbert space is set to be D = 100. For short time 
simulation we set J8t = 0.005 in the Trotter slicing while 
for effective dynamics we set JSt = 0.01 — 0.02 to balance 
the Trotter error and truncation error. Similar procedure 
has recently been used to investigate the decoherence and 
entanglement dynamics induced by spin bath. We hence 
refer to Ref.[27j and the references therein for details of 
simulation procedure. 

From U (t) , one obtains the reduced density matrix p r 
of qubits by tracing out the environment 



p r (t) = Tr bath U(t)p(0)U^(t) 

= ^($ E \U(t)p(p)Ui(t)\$ E . 



(10) 



Here \<&e) is a complete set of state for the spin chain and 
p(0) is the total density matrix at t = 0. Since the initial 
total wavefunction is a product of state, one may consider 
p(0) = \p)(q\ <B> \G)(G\, where \p) (\q)) is an eigenstate to 
the total s z of qubits with the eigenvalue being p (q). At 
time t, the reduced density matrix is given by 



P. 



.(*) - 2<*^l £7 (*)b)|G>(G ! K«l^ t (<)|^). (11) 



For Ising-like coupling, h will be taken to be x so that 
one can replace a z by ±. Therefore, one has 

U{t)\±) = e - l ( t -*«)(tfo±(-l) N eS*)... e - J (t 2 -t 1 )(// T^) 

= U ± (t)\±). (12) 
Hence one can replace U{t)\p) by U p (t)\p). We find 



Pr(t) = \p)(q\ 



Y,($E\U p (t)\G)(G\UtmL 



L E 



(13) 



Since U p (i) no longer acts on \p), one can switch the order 
of (<f> E \U p (t)\G) and (G\U$(t)\$ E )- Using the complete- 
ness of |$e), one obtains 



p r {t) = \p)(q\(G\UUt)U p {t)\G). 



(14) 



Hence the matrix element of the reduced density matrix 
is given by 



p^{t) = {G\UUt)U p {t)\G). 



(15) 



It is clear that the effectiveness of DD control for the 
Ising-like coupling is determined by U%(t)U p (t). 

On the other hand, if the coupling between the qubit 
and the spin bath is Heisenberg-like, \p) is no longer 
an eigenstate to U(t). The reduction from Eq.([TU)) to 
Eq. (fT4l is generally not possible except for the diagonal 



elements^. Therefore, one resorts to Eg. ([10")) to calcu- 
late the reduced density matrix. In this case, the matrix 
element of the reduced density matrix is given by 



P?(t) = Tr hath (p\U(t)p(0)U\t)\q) 
= Tr hath [p(0)UHt)\q)(p\U(t)] 



(16) 



In general, \q) (p\ does not commute with U (t). Hence the 
effectiveness of bang-bang control for the Heisenberg-likc 
coupling is determined by the evolution operator U (t) . 
Note that for longitudinal component when p = q = ± 
for a single qubit, Yang and Lh>i& noticed that for 
the UDD-N pulse sequence applied at a single qubit, 
U(T N ) = ex V [iH eff T N + 0{T% +1 )] and H eff com- 
mutes with |+)(+| and |— )(— |. As a result, the lin- 
ear term in T/v gets canceled and thus the magnetiza- 
tion (a z (T)) can be controlled to 1 + 0(T^ +1 ). Appar- 
ently, the same cancelation does not happen for the off- 
diagonal matrix elements where H e f f does not commute 
with Therefore, to find an effective sequence for 

both dephasing and longitudinal relaxation, one needs to 
directly control U it) to the required order. 



III. SHORT TIME AND LONG TIME 
DYNAMICS OF SINGLE QUBIT 
DECOHERENCE 

In this section, we examine dynamics of a single qubit 
coupled to the spin bath with and without decoupling 
pulses. For short time dynamics of a qubit under DD 
pulses within a cycle of T, we construct pulse sequences 
for suppressing lower orders of U(T) up to 0(T 4 ). Dif- 
ferent strategies for inserting computing pulses are com- 
pared. By using t-DMRG, we demonstrate that comput- 
ing after decoupling performs the best. Therefore, we 
extend DD over a cycle of T to the periodic scheme in 
which computing pulses are inserted at nT. In partic- 
ular, we shall compare long time dynamics of different 
sequences at nT to determine the optimized sequence. 



A. Ising-like coupling 

We start by examining the case when the qubit-bath 
coupling is Ising-like. In this case, there is no longitudi- 
nal relaxation since [s z ,H] = 0. The pure dephasing is 
characterized by the Loschmidt echo L(t) =| p+~(t) | 2 , 
where p\ is defined by Eq. (fT5|) . In the absence of decou- 
pling pulses, it is known that the Loschmidt echo decays 
as L(t) = e~ Q * for short times^. Hence a characterizes 
the short time decoherence of a single qubit in spin bath. 
When the spin bath is modeled by a XY model(i.e., A = 
0) and Hi n t = cs z S z , a can be exactly calculated and will 
be served as a checking point of our t-DMRG numerical 
code. To find a, we first note that after the Wigner- 
Jordan transformation, the Hamiltonian is quadratic and 
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is given by H = J J2 n <4 c n+i + h.c. + es z c\ci . As a re- 
sult, L(t) can be expressed as22 



L(t) =| det[l +r(e 



itH - — itlT* 



-1)] 



(17) 



Here if ± are the matrices corresponding to Hq with s z = 
±1. r is a 2N x 2N matrix with N is the length of the 
spin chain and its element is given by ry = (c\cj). By 
using the identity for the operator A 



det 



l + A 



a Tr ln(l+i) 



TrA-TrA 2 /2+- 



(18) 



and expanding e ltH = 1 + itH ± - t 2 (H ± ) 2 /2 • • • , one 
find a = e 2 . 

On the other hand, in the presence of N decoupling 
pulses in a period T/v, it has been shown that the UDD-N 
pulse sequence suppresses^ (Tn)U+(Tn) = 1 + 
0(T^ +1 ) and the effectiveness increases as N increases. 
As a check, in Figj^l we show our numerical simulations of 
L(T) versus T using t-DMRG. For the free decay without 
decoupling pulses, fitting within < Jt < 0.1, we find 
that a = 0.0224 for J = 1 and e = -0.15. This is 
in agreement with the analytic result 0.0225 within the 
error caused by the Trotter time slicing (J6t = 0.005). 
We also observe that higher order UDD pulse sequence 
is more effective as expected. 




FIG. 2: (Color online) Numerical simulations of Loschmidt 
echo L(T) versus T by t-DMRG for Ising-like coupling. Here 
for a given T, all pulses for any pulse sequence are arranged 
within T. The fitted decaying parameter a for the free case 
is 0.0224 , in agreement with the analytic result. It is clear 
that higher order UDD pulse sequences are more effective. 
Note that in this and the following figures, error bars are 
given for each data point with figure legends being positioned 
at the intersection point of the curve with the error bar. In 
the current figure, errors are about the same sizes of figure 
legends. 



B. Short-time behavior of Heisenberg-like coupling 

When the coupling of the qubit to the spin bath is 
Heisenberg-like, s z is no longer a good quantum number. 
To suppress both longitudinal and transverse relaxations, 



we consider N 7r-pulses with alternating directions along 
x and z axes in the period of T/v . The evolution operator 
then becomes 

• • • a.e^^^-^V.e-^ * 1 [l + 0(Ne Jt 2 )] . 

(19) 



Since for any operator O, one has 

o a O{d)<j a = 0(a a aa o 



(20) 



By inserting appropriate identities, er 2 = 1, one can move 
all the spin operators to the left and obtains 

UW( t ) = {- l ) N a x a z ---a z a x U^ N) [l + 0{NejT 2 p )\ 

= H)V p f/ (Ar) [l + O(iVeJr p 2 )]. (21) 

Here p is an integer and has no contribution in p r . 
Op is a spin operator representing the net operation by 
Vx&z ■ ■ • Vz&x- For example, when N = 3, a p = a z . It is 
clear that a p acts as a parity pulse that compensate fast 
changes due to pulses. To remove its effect, we must add 
an additional pulse a p at T so that U^ N \t) — Uq (t) 
where Uq N ^ is the evolution operator for ideal pulses 
without width and is given by 



U^ N) {t) = e - iH o aN+1 {t-tx) e -iH^ a » (t N -txr-x) x 



. e -iH* 3 (.t 3 ~t 2 ) e -iH* 2 ( t2 - tl ) e -tHt f 1 t x _ 



(22) 



In this expression the time-dependent effective Hamil- 
tonian is defined as 

= H bath + e -[f x (t)o x S* + f y (t)a y Sf + f x {t)*,Sf], 

(23) 

with fp = ±1, /i = x,y,z depending on the sequence. 
Consequently, 

U^ N) (t) = e ~ lHb ^ ht x 

= e - tHbatht G N (t). (24) 

Here f x ,y,z{s) characterizes the sequence history of signs, 
while S^(t) is the operator in the interaction picture: 



St J -{t)=e iHbMht S^e- iHh ^ t 

n=0 



^ j— [Hbath, [Hbath, ' ' 1 [Hbath, 5f] ■ ■ • ]] 



-folds 



E S n tn - 



(25) 



n=0 



It is clear that since H ba th commutes with the qubit spin, 
e -iH bath t j n Eq jjlMj) gets canceled in Eq. dTgl) Hence one 
only needs to suppress G N (T) to the desired order. For 
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this purpose, we re-express G (T) = exp(0 (T)) and 

use the Magnus expansion^ Cl N (T) = Of + 0^ H to 

evaluate G N (T), where by setting A(s) = f x (s)<r x S x (s) + 
fy(s)a y S y (s) + f z (s)a z S z (s), one has 



Since [c^Sf , a v S" ] = 0, all double integrals vanish and 
only single integrals contribute. There are three indepen- 
dent operators in Eq. (|26|) . Consequently, in addition to 
Eqs.(J25]), we require the 1st moment of / M to vanish 



Of (T) 
Of(T) 



A(s)ds 



(26) 



ds 1 / ds 2 [A( Sl ),A(s 2 )]. (27) 
Jo 



If we want to suppress the decoherence to m-th order 
we must suppress Uq N \Tn) to 0(T™), but at the same 
time due to the finite width of pulses we also need to 
ensure O(Nt^) < 0(T"). Using Eq.©, we find the 
minimum of T/v is determined by 



Tjv > max 
> max 



(N + 1)t p ,(NJ 



^,(jVJ 2 -"V B 2 ) 1/m 

rv - " 



2-m Jil/r 



(28) 



We note that when m is small, the first term domi- 
nates. It is thus sensible to define a minimum period 
as T c — Tp/min(oii). Another important observation is 
that due to the finite width of pulses, increasing number 
of pulses also increase TV which leads to stronger deco- 
herence. To find long-time dynamics at t* = nTpf, we 
shall start from m = 1 and focus on small m. For a fixed 
m we find the minimum number of pulses needed, iden- 
tify the optimized value of ctj, and compare the results 
from different m. Indeed, as we shall see, increasing m 
degrades the long-time dynamics in some cases. 

To suppress G N (Tjv) to 0(Tjy ), we keep the mth order 
term in Eqs. (l26t and (l27l) . For m = 1, there are three 
independent operators, (J^S^ (/x = x,y, and z), whose 
coefficients have to vanish. We hence obtain three con- 
straints 



U(t)dt = o. 



(29) 



These constraints can also be obtained from a geometric 
perspective, as pointed out in Refl3ll. Therefore we find 
that the minimum number of pulses is N\ = 4 (including 
the parity pulse). Eqs. (|29[) provide three equations for 
intervals on where % = 1, 2, 3, 4. Together with J2i a i = 1j 
we find that there is only one solution with alternating 
x-z, equally spaced pulse sequence 



ai = 0:2 



a 4 =j. 



(30) 



For m 



2, there are two additional contributions 



in Eqs.(26]) and ([27]) to 0(T£) in G N {T N ). By using 
Eq. ([2^)l . terms with double integrals can be rewritten as 



h = 



dti / di 2 / M (*i)/ v (*2)[5' t (ta),S , '(ti)]. (31) 



*/ M (t)<ft = 0. 



(32) 



By solving Eqs. ([29)) and (|32|) for all possible directions 
of pulses, we find that there are 5 solutions. To continue 
from the case of m — 1, we shall adopt the alternating x- 
z pulse sequence and refer the reader to the Appendix for 
the remaining sequences. Hence, by including the parity 
pulse (y-pulse) , the minimum number of pulses for m = 2 
is N 2 = 7 with intervals for six pulses being 

7- V33 1 
"i = — 77 — = a 7, ot 2 = - = a 6 , 



16 



Q!3 



33-3 
T6 



015, «4 = ^- 



(33) 



Similarly for m = 3 the second moment of has to 
vanish. 



t 2 U(t)dt = 0. 



(34) 



Additionally, there are six mixed moments which have to 
vanish too 



dti / dt 2 (h - t^if^Mh) + Ut 2 )f v (h)\ = 0. 



(35) 



Therefore, in general, the minimum number of pulses for 
m = 3 is N3 — 15. However, for spin- 1/2 qubits, since 
cr 2 = 1 and (5**) 2 = ft 2 /4, terms with /j, = v in Eq.([35|l 
have no contribution to dynamics of qubits. Therefore, 
the minimum number of pulses for a spin-1/2 qubit is 
A^3 = 12. The equations of ctk imposed by Eq. ([33|) is 
generally very complicated. Therefor, one has to resort 
to numerical methods to obtain solutions. By solving 
Eqs.JH, (J32J), (J34]) and ([35]) for alternating x-z pulse 
sequence, we find the numerical values for intervals are 



ai ^ 0.0171, a 2 
a 5 = 0.1184, a 6 
a g = 0.0823 aio 
ai3 = 0.0318. 



: 0.0468, a 3 = 0.0658 a 4 = 0.1013, 
: 0.1006 a 7 S 0.1195, a 8 ^ 0.1049, 
0.1025 an 0.0647 a i2 = 0.0439, 

(36) 



This is the optimized sequences for suppressing U(T) 
to 0(T 4 ). As indicated in the b cginning, due to finite 
width of pulses, m-th order (or 0(T' m )) sequence is not 
necessarily more effective than the (m — l)-th order (or 
0(T m_1 )) sequence. Hence we shall stop at 0(T A ) and 
compare the performance of different orders at nT. 

In Fig. 3(a), we first show numerical results of L[t) 
for < t < T based on t-DMRG for different sequences 



7 



(a) 
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-0.005 
-0.01 

-0.015 
-0.02 









-0.01 




\ * — * OfT 3 ) sequence 








\ * — * 0(T ) sequence 




P 




\ — QDD-2 sequence 




-0.02 




\ ► — ► No pulse applied 




Oil 






\ • • 0(T") sequence 

\ , 1 








• — ♦ OfT" ) sequence using the ground state of bath 
G-O 0(T" ) sequence using the first excited state of bath 




FIG. 4: (Color online) Comparison of Loschmit echo by using 
different initial states. Note that errors are smaller than sizes 
of figure legends. 



FIG. 3: (Color online) (a) Loschmidt echo L(T) versus T 
of different orders simulated by t-DMRG for Heisenberg-like 
coupling with e = —0.3. Note that errors are about the same 
sizes of figure legends. Clearly, L(T) gets improved for in- 
creasing number of orders, (b) L(T C ) versus Jt p . It is seen 
for as r p decreases, 0(T 3 ) sequence starts to win over 0(T 2 ) 
sequence. 

with a fixed T. This corresponds to the ideal pulse sce- 
nario in which the pulse width is neglected. Here we do 
observe that higher order sequences are more effective. 
To take the finite width into consideration, in Fig. 3(b) 
we show numerical results of L(T C ) versus Jt p for se- 
quence of different order. One should keep in mind that 
T c depended on t p . One can see clearly in Fig. 3(b) 
that as t p decreases, 0(T 3 ) sequence starts to win over 
0(T 2 ) sequence. We note that for a much smaller t p the 
0(T 4 ) sequence will out-perform 0(T 2 ) and 0(T 3 ) se- 
quences (not shown in the figure). As a check on if our 
results depend on the initial state of the spin chain, in 
Fig. 4, we compare the Loschmit echo versus T by using 
|$(0)) = \g) ® \bath) with the state of bath \bath) being 
the ground state or the first excited state of the quantum 
spin chain. One can see that their difference is within the 
error bar. In addition, we can also form entangled state 
of the qubit and the bath to check the performance of 
our sequences. Since after the first cycle T, qubit and 
bath is entangled. Therefore, the performance at later 
times is an indirect check and this will be done in the 
next section. 

In Fig. 5, we compare longitudinal relaxation under 
UDD-3 and our optimized sequences. It is clear that the 
optimized sequences outperforms the UDD-3 sequence. 



C. Long time dynamics of Heisenberg-like coupling 

Given the dynamics of dynamical decoupling, one 
needs a strategy for inserting the computing pulses. For 
this purpose, we compare three possible ways to insert 




FIG. 5: (Color online) Comparison of longitudinal relaxation 
of different orders with UDD-3 and QDD-2 sequences. Here 
e = —0.3 and errors are about the same sizes of figure legends. 



a computing pulse that rotates the qubit by an angle 9: 
(i) using a constant B over T in Eq.Q (ii) applying a 
short pulse to rotate the qubit by 9 at some moment t 
with < t < T (iii) applying a short pulse at T after all 
decoupling pulses. From the construction of 0(T N ) se- 
quences, since H sys in Eq.([T]) can be combined with Hi nt 
with Si being replaced by Si+ B, one expects that com- 
puting pulses inserted in < t < T will be suppressed. 
In Fig. 6, we show that the distance of p r to the de- 
sired reduced density matrix versus 9 for three different 
schemes. Indeed, the results clearly show that DD pulses 
suppress the computing as well if one adopts the scheme 
of computing while decoupling. Hence we shall adopt the 
scheme for inserting computing pulses after decoupling 
pulses. 

In the following, we shall adopt the scheme that com- 
puting pulses are inserted at nT while the same DD 
pulses are applied between computing pulses. We would 
like to demonstrate that in this scheme, there will be 
cross-overs of the relative effectiveness for different se- 
quences. Based on the above optimized sequences within 
a cycle of T c , the performance at nT c can be deduced. 
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0(T") sequence - computing by uniform B over T 

2 

B B 0(T ) sequence - computing at T 

— 0(T 3 ) sequence - computing by uniform B over T 

* * O(T) sequence - computing at T 

^\p- m - p (8)1 

2 

A A 0(T") sequence - computing at 0.35T 
0(T 3 ) sequence - computing at 0.7 T 




0.2 



0.4 



0.6 



o.s 



e/jt 



FIG. 6: (Color online) Comparison of three differen schemes 
for inserting a rotation 9. Note that errors are smaller than 
sizes of figure legends and B = 25 Tesla so that 0(T 2 ) se- 
quence outperforms 0(T 3 ) sequence. Here initially the spin 
qubit is in the state |+) and p(8) is the corresponding density 
matrix by applying a rotation of 9 along x-axis on |+) di- 
rectly. D = \p r — p(6)\ is the distance between the simulated 
reduced matrix and p(0)\ with \A\ = TrV A^A 



Since U(nT c ) = U(T c ) n , we have 

U^ N \nT c ) = {(-iyuW [1 + 0{NeJrl)]} n 

« {{-if [1 + 0{e,J m - 1 T m )]} n 

« {H) p }" [1 + O^eJ 1 ™-^" 1 )] . (37) 

Consequently U(nT c ) is suppressed to 0(nT m ). To quan- 
tify the performance at nT c , it is useful to define the 
number of quantum manipulation by defining N op as 



1 



L(N op T N ) = -L(0) 



(38) 



Intuitively N op represents the maximal attainable num- 
ber of quantum operation before the qubit becomes de- 
coherent. The exact number of N op can be found by 
t-DMRG. In Fig. 7(a) we show the effective dynamics 
for lower orders optimized sequences. It is clear that de- 
coupling pulses do suppress the decoherence of a single 
qubit. In Fig. 7(b), we plot N op versus B for lower order 
sequences. We find that m = 1 sequence is the most opti- 
mized sequence for low magnetic fields. This is consistent 
with estimation in Fig. 3. 



IV. SHORT TIME AND LONG TIME 
DYNAMICS OF TWO QUBITS 

In this section, we extend the construction of pulse se- 
quence to multi-qubits. In particular, we shall examine 
the validity of our construction for two qubits using t- 
DMRG. We start by considering n qubits denoted by a 1 
with i = 1,2, • • • ,n. Following the derivation of Eqs. (l2~Tj) 




(a) 



No pulse applied 
3 0(T") sequence 
a 0(T 3 ) sequence 
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Jt* 



o 




(h) O-O 0(T ) sequence 
O— &0(T ) sequence 
A— ^ 0(T 4 ) sequence 



aaaaaAaaaaaaa aa Aaaa a -^ 



B(T) 



FIG. 7: (Color online) (a) Effective dynamics of L(t*) for 
lower orders optimized sequences shows better coherence than 
that of the original dynamics. Here JSt — 0.01, e = —0.3, 
and t* = t/T c {B) with B = 20T. (b)N op for different orders 
versus B fields. In low fields, 0(T 2 ) sequence is the optimized. 
For modest high fields, 0(T 3 ) sequence starts to win over 
0(T 2 ) sequence. 



and (f2"4"|) . one can move all the spin operators due to the 
decoupling pulses to the left and obtain the correspond- 
ing G N (t). The decoupling sequence generally introduces 
different history of sign characterized by /* (t) . We find 

G N (t) = re- l i^ dsk(s ^ with R(s) is given by 

R(s) = £ [fi(*HS?(s) + P v {s)alSf{s) + fi(s)aiSf(s)] 

If one assumes that different sequence is applied to dif- 
ferent qubit, it is clear that in addition to constraints 
set by Eqs.flMl), (02), fl33} and (JHSJ) , there will be extra 
constraints due to cross product of different qubit spins. 
Therefore, it is clear that minimum number of pulses can 
be achieved by setting all the pulse sequence the same: 
/*(£) = ffj,(t). In this case, we find 

R(s) = /x(s)Qx(s) + fy(s)Q y (s) + f z (s)Q z (s), (40) 

where Q^(s) = J2i a li^i( s )- ^ n comparison to the case 
of a single qubit, here replaces the role of a^S^. The 
commutators of are given 

[Q»(s),Q4s')] = Y. a >iW( s )' S j( s ')]- ( 41 ) 

According to Eq. ([23|) . the 0(s n ) of 5f (s) is a n-fold com- 
mutator of Hbath and Since Hb at h only contains 
couplings between nearest neighboring S'f, to 0(s n ), 
S^(s) contains spin operators up to S°± n . Therefore, to 
0(s n s' m ) (i.e., to 0(T"+ m+1 )), we find [£?(«), S%{s')} = 
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if \i — j\ > n + m. 
become 

[Q»(s),Qv(s')] 



Consequently, commutators of Q p 



E 



ay v [S?(*),Si{*% 



(42) 



It is then clear that for spin-1/2 qubits, because — 
1, only commutators with p, ^ v contribute dynamics. 
Furthermore, coefficients that are associated with these 
commutators are exactly the same as those for a single 
qubit. Hence Eqs.flM]), ([52"]). ([55)1, and are also the 
constraints for two qubits to suppress U(T) to 0(T 4 ). 
In other words, both the entanglement of two qubit and 
decoherence of a single qubit can be optimized by the 
same sequence. 




FIG. 8: (Color online) (a) Concurrence C(T) versus T of dif- 
ferent orders simulated by t-DMRG for Heisenberg-like cou- 
pling with e — —0.3. Clearly, C(T) gets improved for in- 
creasing number of orders, (b) C(T C ) versus Jt p . It is clear 
that as t p decreases, 0(T 3 ) sequence starts to win over 0(T 2 ) 
sequence. 

To check the validity of the above conclusion, we ex- 
amine the entanglement of two qubits. To character- 
ize the entanglement, we shall use concurrence as the 
measurement of entanglement^ For a given reduced 
density matrix p{t), the concurrence is defined as C = 
max{Ai — A2 — A3 — A4,0}, where Ai > A2 > A3 > A4 
are the square roots of the eigenvalues of the operator 
p{a v <X> <J v )p*{a v ® <7 V ) and p* is the complex conjuga- 
tion of p. In Fig. 8(a), we show the concurrence C(T) 
calculated by t-DMRG for various sequences. It is seen 
that the sequence in Eq. (I36[) indeed outperforms other se- 
quence. In Fig. 8(b) we show numerical results of C(T C ) 
versus Jt p for different orders. It is clearly seen that as 
t p decreases, 0(T 3 ) sequence starts to win over 0(T 2 ) 
sequence. In Fig. 9(a), we show the effective dynamics 
for the concurrence at t* = nT . In comparison to the 
case without decoupling pulses applied, it is clear that 
decoupling pulses do improve the entanglement of two 
qubits. Fig. 9(b) shows that except for 0(T 4 ) sequence 



with qubits at (39,40), decoupling pulses also suppress 
the generation of entanglement. One of the reasons be- 
hind the non-suppression of the entanglement generation 
for the 0(T 4 ) sequence is due to the large T c required 
by finite t p . However, as indicated by Eq. (l42|) . the dis- 
tance between qubits is also an important factor. As a 
comparison, C(t*) of the 0(T 4 ) sequence for two-qubits 
located at different distances, (39,40) versus (30,50), are 
calculated. It is seen that entanglement generation is 
suppressed only for qubits located at (30,50), in agree- 
ment with conclusions based on Eq. (|4"2")) . In Fig. 9(c), we 
plot N op for different orders versus B fields. We find 
that at modest magnetic field, the lowest order, to = 1, 
is the most optimized sequence for preserving the entan- 
glement. 
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\ • — 0(T 3 ) sequence 
V > — c No pulse applied 
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FIG. 9: (Color online) (a)Effective dynamics of the concur- 
rency. Here JSt = 0.01, e = -0.3, and t* = t/T c (B) with 
B = 20T. (b)Effective dynamics for the generation of entan- 
glement. Here JSt = 0.02, e = -0.9, and t* = t/T c (B) with 
B = 12T,20T,30T for 0{T 2 ),0(T 3 ), and 0(T 4 ) sequences 
respectively. Note that errors are smaller than sizes of figure 
legends. It is clear that entanglement generation is suppressed 
for qubits that are far apart. (c)N op obtained from C(t*) for 
different orders versus B fields. In low fields, 0(T 2 ) sequence 
is the optimized. For modest high fields, 0(T 3 ) sequence 
starts to win over 0(T 2 ) sequence. 



V. SUMMARY AND OUTLOOK 

In summary, feasibility of decoupling pulses that pre- 
serve the coherence and entanglement of spin qubits with 
general couplings to a quantum spin chain are examined. 
It is shown that error induced by pulse- width accumula- 
tion restricts the consideration to sequences with number 
of pulses being less than 10 within each cycle. Within 
each cycle T, exact nontrivial pulse sequences are con- 
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structed to suppress the qubit-bath coupling to 0(T N+1 ) 
progressively with number of pulses being 4, 7, 12 for 
N = 1,2,3. It is demonstrated that computing after 
decoupling has the best performance. Therefore, the per- 
formance is examined by counting number of computing 
pulses that could be applied periodically at the end of 
each DD cycle. Based on t-DMRG, our numerical re- 
sults show that for modest magnetic fields (10-40 Tesla) 
available in laboratories, the overall performance is opti- 
mized when number of pulses in each cycle is 4 or 7 with 
pulse directions be alternating between x and z. 

While so far our numerical results are obtained by us- 
ing either a single qubit or two qubits as demonstrations, 
our results also provide insights for preserving coherence 
and entanglement of multi-qubits. In fact, according to 
our analysis, in principle the same sequences we obtained 
in this work can also dynamically decouple multi-qubits 
from the environment in low magnetic fields. For high 
magnetic fields, to obtain better coherence and entangle- 
ment, one needs to suppress higher order terms. In this 
case, however, one still needs to suppress lower orders. 
Therefore, our results will still serve as a useful starting 
point for sequences for higher magnetic fields. 
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Appendix A: 0(T 3 ) Optimal Sequences 

In this appendix, we explicitly construct all possible se- 
quences that suppress U(T) to 0(T 3 ). For this purpose, 
we first note that operators appear in U (T) in the order of 
O(T) and 0{T 2 ) are a x S x , a y Sy, a z S z , [H , [H ,S x }}a x , 
[Ho,[Ho,S v ]]a y , and [Hq,[Hq, S z ]]a z . Requiring coeffi- 
cients of these operators to vanish yields Eqs. ([29"l) and 
([32]) . Eqs.d29]) and (02) can be solved by using Math- 
ematica to exhaust all possible configurations of pulse 
directions. We find that in addition to Eq. ([33)l . the fol- 
lowing generic sequences are also solutions 
Sequence 1: pulse direction xzxxzx 
a.% = ot-2 = 013 = ol§ = ot§ = ai = 1/8, Qf4 = 1/4 
Sequence 2: pulse direction xzxxyx 
at = 0.104715, a 2 = 0.145282, a 3 = 1/8, a 4 = 1/4 
a 5 = 1/8, a 6 = 0.145282, a 7 = 0.104715 
Sequence 3: pulse direction xzxzyz 
ai = 0.0785, a 2 = 0.1396, a 3 = 0.1596, a 4 = 1/4 
a 5 = 0.1715, a 6 = 0.0931, a 7 = 0.1104 
Sequence 4: pulse direction xzxyzy 
ai = 1/8, a 2 = 0.095491, a 3 = 0.1545, a 4 = 1/4 
a 5 = 0.1545, a 6 = 0.095491, a 7 = 1/8. 
In addition to the above sequences, equivalent sequences 
can be formed by applying cycle permutations on x, y 
and z. 
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